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' Abstract. We present a combined analysis of the optical spectral variabil- 

^sj" I ity for two s amples of QSOs. 42 objects at z < 0.4 monitored at the Wise 

Observatory ijGiveon et al. llT999|) . plus 59 objects up to z 3 in the field of 
the MaKel lanic Clouds, detec ted and/or monitored within the MACHO Project 
^ I database l|Geha et al. 1120031) . Our analysis shows some increase of the observed 

]^ . spectral variability as a function of redshift, with a large scatter. These data are 

compared with a model based on the addition of flares of different temperatures 
to a stationary quasar SED, taking into account also the intrinsic scatter of the 
SEDs. 
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Introduction 

O 

^ I Variability of the spectral energy distribution (SED) of Active Galactic Nuclei 

c/3 . (AGNs) is a powerful tool to investigate the role of the main emission processes, 

and the origin of their variations. Many variability mechanis ms have been pro- 
posed in the past, including su pernovae explosions (Aretxaga et al. iFlOOT ^ in- 
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stabilities in the accretion disk (iKawaguchi et al. 119981 ) ■ and gravitational lens- 
ing due to intervening matter ([Hawkins Ill983f ). Since it is likely that different 



\ physical phenomena are causing variability in different bands and timescales, 

multifrequency analyses will be ultimately needed to ob tain a complete de- 
scrip tion. However, we have shown in a previous paper ( Trevese Sz Vagnetti I 



that even a two optical band analysis, once performed on a statisti- 
cal sample, provides valuable constraints on the origin of variability. In the 
optical band, AGNs most commonly become bluer, i.e. their spectrum be- 
comes harder, when brighter. This has been shown t o occur for some in- 
dividual AGNs (ICutri et al I Il98,4 lEdelson et aT~l 119901 : iKinnev et afl Il99ll : 
iPalt ani Sz Courvoisier lll994l T! and for a few complete samples: P G quasars at 
z < 0.4 monitored at the Wise Observatory ( Giveon et al. l[l999l . hereafter re- 



ferred as WO quasars), f aint quasars in the SA 57 with two-epoch informa- 
tion (iTrevese et ah Il2n0lh . and quasars from the MACHO Project database 
( Geha et al. Il200,l lereafter referred as MP qua sars), which are here studied 



toget her with the WO quasars (see also IVagnet ti &: Trevese I l2003l ). Our anal- 
ysis ( Trevese &: Vagnetti Il2003 ) of the two-color WO light curves confirms the 



hardening in the bright phase of individual quasars. This in turn accounts 
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for two other effects: (i) average variability of quasar samples is greater at 
higher rest-frame freq uencies (Di Clemente et al. 1996), a results co n firmed by 
recent results of IVand en Be rk et al. I (|2004l '). Ide Vries et al. I (|2003l ). Ide Vries I 
(200l), and (ii) - at a given observed frequency - variability is co rrespondingly 
greater at higher redshifts ((Giallongo. Trevese fc Vagnetti "l99l'). To quantify 
the amount of spectral variability, we iTrevese &: Vagne tti 2002) introduced the 
spectral variability parameter (3 = Aa/Alog/jy {f^, being the specific flux and 
a = dlog fu/dlogiy the spectral slope) and the use of the a-f3 plane to derive 
constraints on the variability mechanisms. In particular, it was shown that the 
variation of the spectral slope implied by changes of the accretion rate is not 
able to explain /3 values deduced from the observations, wh ile "hot spots" on 
the accretion disk, possibly caused by instability phenomena (jKawaguchi et al~l 
can easily account for the observed spectral variability. 



Samples 

The present an alysis is based on the sets of light curves for the Wise Observatory 
(WO) sample (iGiveon et al. Ill999h and for the MACHO Project (MP) sample 
( Geha et al. 11200.-^ . l)oth available in electronic form. The WO sample consists 
of 42 PG quasars selected to be nearby, i.e. z < 0.4, and bright, i.e. i? < 16 mag, 
monitored in the Johnson-Cousins B and R bands with the 1 m Wise Observatory 
telescope for a total duration of 7 years and a median observing interval of 39 
days. The MP sample includes 47 quasars discovered behind the Magellanic 
Clouds on the basis of their variability, and 12 additional quasars previously 
known in the same area. All the 59 MP quasars were monitored for 7.5 years 
with Mount Stromlo Observatory's 1.27 m telescope, with average sampling 
times of 2-10 days. Magnitudes in the special MACHO red and blue passbands 
were then transformed and made available as standard V and R magnitudes. 
The MP quasars cover the redshift interval 0.2 < z < 2.8 and have apparent 
mean magnitudes 16.6 < V < 20.1. The combined analysis of the two samples 
allows the investigation of the spectral variability as a function of z. 



Single Band Variability 

Analysis of the dependence of variability on redshift and luminosity on single 
magnitude-limited samples is affected by the strong L-z correlation intrinsic to 
the samples themselves, therefore it is preferab le to combine more samples to 
cover a larger portion of the L-z plane (see e.g. iGiallongo. Trevese &: Vagnetti I 
[l993). Thus, we combine the MP sample, with data in V and R for 59 quasars, 
and the WO sample, which h as variability information in B and R for 42 PG 
quasars at z ^ 0.4. Following IGiallongo. Trevese fc Vagnetti I (|l99ll ). we mea- 



sure variability in a fixed bin of the rest-frame Structure Function, 75-150 days. 
For the R band, which is common to the two samples, we find no correlation 
of variability with redshift for the MP sample alone, but we find an increase 
of variability with redshift for the combined samples, with correlation coeffi- 
cient r = 0.19 and probability P{> r) = 0.05, and a more significant result, 
r = 0.44, P(> r) = 3 • 10~^, if we consider only objects of the two samples in 
a limited range of absolute magnitudes (— 26 < Mr < —22). We find also that 
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variability increases with M/j, with r = 0.25, P{> r) = 0.01 for the combined 
samples, and with r = 0.30, P{> r) = 0.02 for t he MP sample alone. These re- 
sults c onfirm t he previous trends foun d, e.g., bv lGiallongo. Trevese h Vaenetti I 
((199 ih and by ICristiani et al. I (Il996t). Variability increases also as a function 
of rest-frame frequency (|Di Cl emente et al. "1996'j: results for the MP and WO 
samples have been shown bv IVagnetti &: Trevese . (,2003 ) . 




Figure 1. Left: spectral slope a vs log z for the two quasar samples; crosses: 
WO sample; circles: MP sample. The continuo us curve represents t he predic- 
tion for the average quasar SED determined bv lElvis et al. I l|1994j) . assumed 
independent of z in shape and plotted as a function of the redshift corre- 
sponding to Vrest (which is labeled in the upper part of the diagram) for the 
observed V band. Dashed curves correspond to the 90% percentile SEDs pub- 
lished by the same authors. Right: The spectral variability parameter [3 vs 
log z for the two sa mples. The curves correspond to model predictions for 
the average SED bv lElvis et al. I l)1994|) . perturbed with blackbody flares of 
temperatures T = 2 ■ 10^ K (upper curve), T — & ■ 10^ K (middle curve), 
T = 3 • 10^ K (lower curve). 



Spectral Variability 

As the MP sample includes objects in a large range of redshift, up to z ^ 3, 
the spectral slope a measured at a given observed frequency is in fact the slope 
of the SED at different emission frequencies, increasing with z. Therefore we 
plot, in Figure 1 {leff) a vs z for the two samples, comp ared with a prediction 
for an average SED, as determined bv lElvis et Xl (11994 ). To take into account 



the large scatter of the quasar SEDs, we a lso show the predictions for the 90% 
percentiles reported bv lElvis et al. I (11994 ). The comparison is consistent with 



the assumption of an average SED approximately independent of z in shape, 
although some larger scatter in the spectral slopes appears in the low-z data. 
To perform a similar comparison for the spectral variability parameter /?, as a 
fu nction of z, we requi r e a m odel for the variability of the SED; we adopt, as 
m evese fc Vagnetti l (l2002l ). the addition of a blackbody flare to the average 



SED bv lElvis et al. I (119^ In Fi gure 1 (right), showing (3 vs z for the two 
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samples, some increase of /? with z is apparent, but there is a large scatter. 
The model predictions, for flares of different temperatures added to the average 
quasar SED, account for at least part of the scatter. 



Conclusions 

The study of spectral variability of quasars is a powerful tool to test the validity 
of different emission and variability models against time-dependent data of sta- 
tistical quasar samples. For such purpose, it is necessar y to obtain var iability 
information in at least 2 bands. The samples of quasars bv lGiveon et al. 1 ( 1999), 
with light curves for 42 quasars in B and R, and that bv iGeha et al. I ( 2003) , 
with information in V and R for 59 quasars, are both appropriate for this pur- 
pose. Moreover, the latter sample allows to extend investigation of the spectral 
variability to high z. Since the MP sample is selected on the basis of variability, 
our analysis could be biased by possible systematic differences in the variability 
properties of the two samples. However, since the variability threshold adopted 
for detection in the MP sample is low (0.05 mag), t he incompleteness res pect, 
e.g., to color-selected samples becomes negligible (see iTrevese et al. 1119891) . In- 
vestigation of spectral variability as a function of redshift corresponds to the 
study of the perturbations of the SED at different rest-frame frequencies. Large 
samples with multi-band variability information are ultimately needed to char- 
acterize the dependence of the spectral variability parameter f3 on both Vrest 
and z. Progr ess in this direction will be achieved by a new analysis of the type 
presented bv lVanden Berk et al. I ( 20041 ) . if repeated observations will provide a 
better time sampling. 
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